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ABSTRACT 

In  this  paper,  we  employ  a  homogenized  energy  model  (HEM)  for  shape  memory  alloy  (SMA)  bending  actuators. 
Additionally,  we  utilize  a  Bayesian  method  for  quantifying  parameter  uncertainty.  The  system  consists  of  a  SMA 
wire  attached  to  a  flexible  beam.  As  the  actuator  is  heated,  the  beam  bends,  providing  endoscopic  motion.  The 
model  parameters  are  fit  to  experimental  data  using  an  ordinary  least-squares  approach.  The  uncertainty  in 
the  fit  model  parameters  is  then  quantified  using  Markov  Chain  Monte  Carlo  (MCMC)  methods.  The  MCMC 
algorithm  provides  bounds  on  the  parameters,  which  will  ultimately  be  used  in  robust  control  algorithms.  One 
purpose  of  the  paper  is  to  test  the  feasibility  of  the  Random  Walk  Metropolis  algorithm,  the  MCMC  method 
used  here. 
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1.  INTRODUCTION 

Shape  memory  alloys  (SMA)  are  unique  actuators  that  are  gaining  increased  utilization  in  prototypes  such  as 
robotic  catheters,1,2  robotic  hands,3  jet  chevrons,4,5  and  underwater  vessels.6  SMA  actuators  are  capable  of 
recovering  large  strains  (approximately  5%)  upon  heating.  However,  the  control  and  optimization  of  SMA  devices 
is  complicated  by  the  material’s  nonlinear,  hysteretic  dependence  on  stress  and  temperature. 

In  this  paper  we  employ  a  computationally  efficient  model  for  a  SMA  bending  actuator.  Furthermore, 
we  utilize  Markov  Chain  Monte  Carlo  (MCMC)  methods  to  quantify  parameter  uncertainty  in  a  systematic 
manner.  MCMC  methods  are  based  on  Bayes’  rule,  which  relates  the  posterior  density  of  model  parameters  to 
a  prior  density  and  the  likelihood  of  those  parameters.7,8  MCMC  methods  avoid  the  difficulty  of  calculating  the 
likelihoods  directly  (and  integrating  over  the  parameter  space)  by  creating  a  Markov  Chain  where  the  stationary 
distribution  is  the  posterior  distribution.  A  proposal  function  and  an  accept-reject  rule  determine  the  random 
walk  of  the  Markov  Chain.9  The  main  goal  of  this  paper  is  to  test  the  feasibility  of  the  approach  for  the  SMA 
bending  actuator  under  consideration. 

An  algorithmic  approach  to  quantifying  parametric  uncertainty  is  useful  for  numerous  reasons.  The  results 
can  be  used  to  produce  bounds  on  the  model  output  for  comparison  to  experimental  data.  An  additional  use  for 
SMA  actuators  (in  particular  the  bending  actuator  under  consideration)  is  to  quantify  parameter  uncertainty  for 
control  algorithms.  Robust  control  methods  require  some  measure  of  model  uncertainty  in  order  to  guarantee 
stability  or  optimality.  For  example,  sliding  mode  control  (or  variable  structure  control)  requires  bounds  on 
model  uncertainty  in  order  to  ensure  the  attractiveness  of  the  sliding  surface.10 

The  remainder  of  this  paper  is  organized  as  follows.  The  model  of  a  flexible  beam  actuated  by  a  single  SMA 
tendon  is  presented  first.  The  homogenized  energy  model  (HEM)  is  used,  which  has  been  applied  to  various  smart 
materials.11  20  The  approach  to  quantify  uncertainty  is  described,  including  optimization  of  model  parameters 
and  the  Markov  Chain  Monte  Carlo  method  used  here.  Finally,  the  results  are  presented. 
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2.  HOMOGENIZED  ENERGY  MODEL  OF  A  SINGLE-TENDON  SHAPE  MEMORY 

ALLOY  BENDING  ACTUATOR 


The  system  consists  of  a  flexible  beam  actuated  by  a  single  shape  memory  alloy  (SMA)  tendon,  as  depicted  in 
Figure  1(a).  The  SMA  tendon  (actuator)  is  held  a  fixed  distance  from  the  neutral  axis  of  the  flexible  beam  using 
rapid-prototyped  collets.  The  SMA  actuator  is  pre-strained  before  it  is  attached  to  the  axially  stiff,  laterally 
compliant  beam.  Therefore,  a  moment  is  created  as  the  SMA  tendon  contracts  due  to  Joule  heating  and  the 
structure  bends  in  the  manner  shown  in  Figure  1(b).  Similar  robotic  systems  include  catheters  and  smart 
inhalers.21,22  A  mesoscopic  free  energy  model  of  a  catheter  was  introduced  in  Veeramani  et  a/.,1  with  further 
work  in  Crews  et  al .2,2S 


2.1  System  Model 

As  shown  in,2,23  the  bending  angle  6{t)  can  be  related  to  the  SMA  tendon  stress  a(t)  by  the  relation 


e(t) 


aAcLcr(t ) 
El 


(i) 


where  a  is  the  actuator  offset  from  the  neutral  axis,  Ac  is  the  cross-sectional  area  of  the  SMA  actuator,  L  is 
its  length,  E  is  the  elastic  (Young’s)  modulus  of  the  flexible  beam,  and  I  is  its  area  moment  of  inertia.  The 
stress-strain  relationship  for  the  bending  system  can  be  simplified  to 


a(t)  = 
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(2) 


where  Ep  is  the  pre-strain  in  the  SMA  actuator.  Equation  (2)  reduces  the  planar  bending  problem  to  a  one¬ 
dimensional  problem  consisting  of  a  SMA  actuator  in  parallel  with  a  linear  spring  with  stiffness 
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The  SMA  strain  e(t)  is  modeled  using  the  homogenized  energy  model.  The  HEM  accounts  for  material 
inhomogeneities  and  interaction  effects  to  accurately  and  efficiently  quantify  the  macroscopic  SMA  behavior. 
A  complete  derivation  of  the  model  and  data-driven  techniques  for  estimating  model  parameters  are  provided 
in.17  In  Crews  et  al .,17  a  single  SMA  actuator  in  constant-temperature  conditions  was  considered.  The  relevant 
relationships  and  equations  are  summarized  here  in  order  to  tailor  the  HEM  to  the  bending  system  under 
consideration. 

The  HEM  quantifies  the  macroscopic  strain 
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Figure  1:  (a)  Flexible  beam  actuated  by  single  SMA  tendon;  (b)  Time-lapse  photograph  of  robotic  system. 


by  assuming  the  relative  stress  a r  and  interaction  stress  07  are  manifestations  of  underlying  densities  vr(cfr) 
and  z/j(ctj).  In  Equation  (4),  T(t)  is  the  temperature  in  the  SMA  wire.  Different  forms  can  be  assumed  for  the 
densities,  subject  to  certain  conditions.11  Here,  we  assume  the  densities 
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vi{ai)  =  7T  E^fc^fe  (<7/)  C2  =  E^fc  (6) 
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are  linear  combinations  of  log-normal  and  normal  functions.  The  kernels  <j)k,m(&R )  and  ^pk{^i)  are  given  by 
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and 
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The  density  coefficients  and  /?*.  and  parameters  C&,  /iRm,  and  bk  can  be  estimated  using  experimental  data.17 
The  mesoscopic  (or  local)  strain  e(a(t)  +  07,  T(£);  ctr)  in  Equation  (4)  is  given  by 

e(oV)  +crr,T(i);crR)  =  xA(t)^E-  +  xM+(t)  ( +  £t)  +xM-(i)  ( “  £t)  , 
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where  x^(t),  and  XM-(t )  are  the  austenite  (A),  martensite  plus  (M+),  and  martensite  minus  (M— ) 

phase  fractions,  respectively,  is  the  austenitic  elastic  modulus,  Em  is  the  martensitic  elastic  modulus,  and 
Sr  is  the  maximum  recoverable  strain. 

The  evolution  of  the  phase  fractions  is  governed  by  the  coupled  differential  equations 


XM+{t)  =  ~p+AXM+(t)  +  pA+XA(t) 

XM-(t )  =  ~P-AXM-(t)  +  pA-XA(t) 

(7) 

and  the  conservation  relation 

Substituting  (8)  into  (7)  yields 

xA{t)  =  1-  xM+{t)  -  xM-(t). 

(8) 

XM+{t) 

=  ~  (p+A  +  PA+)  XM+(t)  -pA+XM-(t)  +pA+ 

=  ~  ( P-A  +PA-)xM-(t )  ~PA-XM+(t)  +PA-. 

(9) 

The  transition  rates  pap  (a,  /3mA,  M+,  M— )  depend  on  the  SMA  stress,  temperature,  and  material  parameters 
and  are  described  in  greater  detail  in  Crews  et  al17  The  temperature  in  the  SMA  actuator  is  quantified  using  a 
first-order,  lumped  capacitance  model, 

t{t)  =  -h{T{t)-T00)+1v\t),  (10) 

where  h  is  the  convective  heat  transfer  coefficient  between  the  actuator  and  ambient  air  at  temperature  X^,  v(t) 
is  the  applied  voltage,  and  7  is  a  parameter  that  accounts  for  the  SMA  resistance  and  specific  heat. 

The  integral  (4)  is  discretized  using  4-point  Gaussian  quadrature  on  20  equal  intervals,  yielding 
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where  and  07.  are  the  quadrature  points  and  Vi  and  Wj  are  the  quadrature  weights.  The  summations  in  (11) 
can  be  converted  to  the  vector-matrix- vector  product 


e(t)  =  vTr  w, 


where 

VT  =  [viuR(aRl),  ■  ■  ■  ,vNivR((TRNi )] 

WT  =  [wii/i((Tix),  ■■■  ,  wNjvi(oiN. )] 

V  =  XA^+XM+  (^+£T)+XM-(S-er). 

The  Ni  x  Nj  matrices  Xa,  Xm+,  and  Xm-  evaluate  the  phase  fractions  at 

[Xa(t,a(t),T(t))]ij  =xa  +  aI:j,T(t)-,aRi)  a  =  A,M+,M~, 

where  xa  are  the  solutions  to  (7)  and  (8). 

The  ODEs  (9)  are  discretized  and  solved  using  an  implicit  Euler  scheme.  For  discretized  time  tk 
stress  ak  =  a (tk),  and  temperature  Tk  =  T(tk)  values,  the  implicit  Euler  discretization  yields 
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Note  that  in  (13),  the  transition  rates  depend  on  crfe,  Tfe,  07,  and  ctr.  The  austentic  phase  fraction 
by  the  conservation  relation  (8).  The  implicit  Euler  discretization  of  the  temperature  yields 
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The  stress  V  is  found  by  substituting  (11)  into  (2)  and  solving  for  the  equilibrium  stress,  which  yields 
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Equation  (12)  can  be  solved  using  Cramer’s  rule,  which  gives 
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(13) 

is  given 


(14) 


The  computational  efficiency  of  the  model  can  be  improved  by  storing  4-D  arrays  Vap  of  the  transition  rates 
evaluated  at  the  quadrature  points  a r.  and  ajj  and  uniformly  distributed  values  of  the  stress  a£  and  temperature 
"T  i  where  crmjn  ft  c  ft  c^max  and  Tmin  ^  T  ft  Tmax.  For  example, 

[P{mA]ij=P+A  (<7t  +  <7Ij,Tm-,<TRi). 

During  implementation,  the  indices  t  and  m  corresponding  to  values  nearest  to  ak  and  Tk  are  found  and  the 
Ni  x  Nj  matrices  are  used  to  calculate  the  matrices  corresponding  to  (14),  C^1,  •  •  •  ,  Ck^x .  These  matrices 
can  be  rapidly  calculated  using  point-wise  operations.  The  matrices  of  the  phase  fractions  can  then  be  updated 
in  one  step  using  the  relations 
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where  point-wise  multiplication  (.x)  and  summation  are  again  used.  In  Equation  (15),  11  is  the  Ni  x  Nj  identity 
matrix. 


3.  PARAMETER  UNCERTAINTY  QUANTIFICATION 

Numerous  techniques  exist  to  estimate  uncertainty  in  model  parameters.  Frequentist  approaches  involve  multiple 
measurements  or  samples  of  experimental  data.  However,  collecting  large  sets  of  data  can  be  expensive  and  time 
consuming.  One  method  that  overcomes  these  limitations  is  bootstrapping.24  26  Here,  we  are  testing  the  Bayesian 
framework  by  using  Markov  Chain  Monte  Carlo  (MCMC)  algorithms.  First,  we  determine  the  optimal  model 
parameters  using  standard  least-squares  approaches.  Then,  we  use  the  optimization  results  to  determine  the 
covariance  matrix  for  the  MCMC’s  proposal  function. 

3.1  Ordinary  Least  Squares  Fit  of  Model  Parameters 

The  optimization  algorithm  minimizes 


m  =  as) 

Z  i= 1 

the  sum  of  squared  error  between  the  experimentally  measured  bending  angle  y(U)  and  the  model  predicted 
bending  angle  y{U\p)  =  6(ti\p)  given  by  Equation  (1).  A  complete  list  of  the  model  parameters  p  is  provided  in 
Table  1.  Here,  we  are  including  all  the  SMA  model  parameters  in  the  optimization  routine.  Alternatively,  one 
could  optimize  the  SMA  model  parameters  (such  as  ak,m  through  V  in  Table  1)  using  tensile-test  data17  and 
then  only  fit  the  heat  transfer  and  bending  model  parameters  using  the  measured  bending  angle.  Including  all 
the  parameters  will  produce  a  lower  sum  of  squared  error  but  neglects  the  coupling  between  parameters. 


3.2  Markov  Chain  Monte  Carlo  Methods 


MCMC  algorithms  are  based  on  Bayes’  rule,  which  relates  the  posterior  density  7r(q\y)  to  a  prior  density  tt pr(q) 
by  the  relation 


n(q\y) 
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(17) 


where  p(y  \  q )  is  the  likelihood  of  observing  y  given  parameters  q.  For  relatively  few  parameters,  the  posterior 
density  can  be  calculated  directly  by  integrating  the  denominator  in  Equation  (17).  However,  in  higher-dimension 
parameter  spaces,  this  approach  is  infeasible. 


MCMC  methods  avoid  the  difficult  integral  in  Equation  (17)  by  creating  a  Markov  Chain  whose  stationary 
distribution  is  the  posterior  density  7r(q\y).  Various  MCMC  algorithms  exist,  including  the  popular  Metropolis- 
Hastings  algorithm7  9,27  and  adaptive  algorithms.28 


Table  1:  Model  parameters  and  estimation  techniques. 


Variable 

Description 

Units 

& k,m 

Relative  stress  density  coefficients 

- 

Pk 

Interaction  stress  density  coefficients 

- 

Ea 

Elastic  modulus  of  austenite 

GPa 

Em 

Elastic  modulus  of  martensite 

GPa 

°L 

Martensite  transition  stress  at  temper¬ 
ature  Tl 

MPa 

Tl 

Lower  transition  temperature 

K 

Actt 

Hysteresis  loop’s  temperature  depen¬ 
dence 

MPa/K 

et 

Maximum  recoverable  strain 

% 

T 

Relaxation  time 

s 

V 

Layer  volume 

m3 

h 

Convection  coefficient 

- 

7 

Heat  transfer  parameter 

- 

£p 

SMA  actuator  pre-strain 

%n 

a 

SMA  actuator  offset  from  the  neutral 
axis 

mm 

L 

Flexible  beam  length 

mm 

El 

Beam  elastic  modulus  and  area  moment 
of  inertia 

N-cm2 

Here,  we  are  using  the  Random  Walk  Metropolis  algorithm  presented  in  Solonen.9  The  algorithm  calculates 
new  parameters  q  new  from  a  proposal  distribution  /(•  |  q  old).  The  proposal  distribution  is  taken  to  be  Gaussian; 
therefore,  we  need  an  estimate  for  the  covariance  matrix  C .  An  initial  set  of  parameters  g*  is  determined  first 
using  standard  least-squares  approaches.  Using  the  assumption  that  the  residuals  are  i.i.d.,  the  covariance  matrix 
is  approximately 

C  =  cr2(JTJ)-\  (18) 

where  J  is  the  Jacobian  of  the  model  error 


e(ti)  =  y{ti)  -y(ti',q*). 

The  variance  in  the  model  error  is  estimated  by 

j2  SS(q*) 

N  —  length(g  )  ’ 


where  the  sum  of  squared  errors  is 


N 

SS(q)  =  ^2  ( y(tt )  -  y(U;  qf 

i= 1 


for  N  experimental  data  points.  After  initializing  the  algorithm  with  the  starting  point  q  *  and  covariance  matrix 
C,  a  new  set  of  parameters  is  proposed  using  the  random  walk  relation 


qnew  =  qM  +  Rz, 


where  R  is  the  Cholesky  decomposition  of  C  and  z  is  a  random  vector  sampled  from  the  standard  normal 
distribution. 


The  new  set  of  parameters  is  accepted  with  probability 


.  (  ir(qnew  |  y)f($old  I 

a  mm  ^  ,  ^old  |  y}f(fnew  |  fold}  J  ■ 


(19) 


Since  the  Gaussian  yields  a  symmetric  proposal  f(qold  \  qnew )  =  f(qnew 
reduces  to 


n(qnew  |  y) 
it (qold  |  y) 


P(y  I  qnew)7rpr(qnew ) 
p(y  I  qold)7Tpr(qold) 


qold ),  the  ratio  in  Equation  (19) 

(20) 


Note  that  in  Equation  (20),  the  constant  JRdp(y  \  q)^Pr(q)dq  cancels  and  we  know  the  posterior  density  up  to 
this  normalization  constant.  Finally,  the  assumption  of  a  uniform  (or  uninformative)  prior 


7 Tpr(q°ld)  =  *Pr(qneW)  =  1, 


yields  the  acceptance  ratio 


Here,  the  likelihood  is  given  by 


a  =  min 


1  p{y  I  <fwew) 


p(y  I  qold) 

p(y  I  q)  =  exp  (-0.5a2  SS(q)2) 


(21) 


with  the  assumption  of  Gaussian  model  errors.  For  the  MCMC,  we  are  using  a  subset  of  the  parameters  listed  in 
Table  1  and  therefore  have  denoted  the  MCMC  parameters  with  q  (versus  the  optimization  parameters  p).  Since 
the  goal  is  to  use  the  uncertainty  quantification  results  in  robust  control  algorithms,  the  heat  transfer  parameters 
h  and  7  are  included,  as  the  heat  transfer  dynamics  has  a  strong  effect  on  the  response  of  the  SMA-actuated 
beam.  Additionally,  the  parameters  related  to  the  bending  model  are  included:  a,  £p,  and  EL  However,  all 
the  parameters  Table  1  can  be  included  in  the  MCMC  without  affecting  the  convergence  of  the  Monte  Carlo 
algorithm,  since  the  convergence  depends  on  the  number  of  iterations  and  not  the  number  of  parameters.  The 
Random  Walk  Metropolis  algorithm  for  the  SMA-actuated  beam  is  summarized  in  Algorithm  1.  During  the 
optimization  and  MCMC  steps,  all  the  parameters  are  scaled  so  that  the  magnitudes  are  on  the  order  of  unity. 


4.  RESULTS 

4.1  Ordinary  Least  Squares  Fit  of  Model  Parameters 

A  prototype  of  a  SMA-actuated  flexible  beam  was  constructed  to  gather  experimental  data.  The  system  consists 
of  a  0.127  mm  diameter  FLEXINOL  SMA  actuator  (Dynalloy,  Inc.  Tustin,  CA)  and  a  0.5  mm  diameter  super¬ 
elastic  Nitinol  beam.  The  bending  angle  is  measured  using  a  trakSTAR  3D  magnetic  tracking  system  (Ascension 
Technology  Corporation,  Burlington,  VT).  Complete  details  of  the  experiment  are  provided  in  Hannen  et  al29 

An  input  voltage  consisting  of  a  sinusoidal  function,  ramp  input,  and  step  input  of  different  magnitudes  is 
used  to  collect  the  experimental  data  shown  in  Figure  2(a).  This  input  captures  the  major  loop  and  various 
minor  loops  for  the  SMA  actuator.  A  comparison  between  the  model  output  using  initial  parameter  estimates 
and  the  experimentally  measured  bending  angle  is  shown  in  Figure  2(b).  A  comparison  between  the  fit  model 
and  measured  bending  angle  is  shown  in  Figure  2(c).  A  complete  list  of  the  initial  parameter  estimates,  bounds, 
and  least-squares  optimal  parameters  is  summarized  in  Table  2.  The  bounds  are  chosen  to  keep  the  parameters 
within  a  reasonable  range  of  the  parameters  found  in  Crews  et  al.17  or  within  reasonable  ranges  based  on 
fabrication  tolerances. 

4.2  Markov  Chain  Monte  Carlo 

The  MCMC  algorithm  is  parallelized  using  8  separate  threads.  Each  thread  is  initialized  with  the  same  starting 
value  (the  optimal  parameters  in  Table  2)  and  run  for  50,  000  iterations  to  ensure  burn-in.  The  parameter  values 
for  each  iteration  are  shown  in  Figure  3.  The  resulting  histograms  are  shown  in  Figure  4. 

The  results  indicate  that  the  optimal  parameters  were  not  obtained  during  the  ordinary  least-squares  mini¬ 
mization.  The  MCMC  eventually  approaches  the  upper  bound  on  7  (5.0)  and  a  convection  coefficient  h  around 


Find  p*  =  minX^i  (y(U)  -  y(U-,p))2  ; 

o'2  =  X)£Li  {y{ti)  -  y(U;P*)f  /  (N  -length(p)); 

ss°ld  =  zZi  (m-y(ti-,p*))2-, 

qold  =  [a* , £p,  ft* , 7* ,  EI*]t  ; 

Calculate  Jacobian  Jq  using  a  finite  difference  method; 

C  =  a2{jTjq)-'  ; 

Calculate  R  so  that  C  =  RT R  (using  Cholesky  decomposition)  ; 

for  i  =  1  to  Ni  do 

Sample  z  =  Zk,  k  =  1, • *  •  ,  length from  iV(0, 1)  : 


qnew  =  qold  +  Rz  ; 

Sample  from  C/[0, 1]  ; 

55-“  =  E,=i  (»(*i)  -  y(U;  q  new,p*))2  ; 
a  =  min {1, exp  (  -  0.5a~2(SSnew  -  SSold ))}; 

if  ua  <  a  then 

_  rfnew  . 

q  q  , 

old  _  T*new . 

q  q  i 


II  7V(0, 1)  is  the  standard  normal  density 
//  U( 0,1)  is  the  uniform  density  on  [0  1] 


gg°ld  _  ggnew. 


else 


_  rfold  . 

q  q  i 


end 


end 


Algorithm  1:  Random  Walk  Metropolis  algorithm  for  the  SMA  bending  actuator. 


Tims  ]s]  Tims  ]$]  Time  ]s] 


(a)  (b)  (c) 

Figure  2:  (a)  Input  voltage  and  bending  angle  comparison  between  model  and  experimental  data  for  (b)  initial 
parameter  estimates  and  (c)  optimal  parameters. 


1.65.  The  root  mean  squared  error  (RMSE)  corresponding  to  these  values  (7  =  5.0  and  h  =  1.65)  is  3.93°, 
compared  to  a  RMSE  of  3.98°  for  the  OLS  parameters  listed  in  Table  2.  The  histograms  of  the  heat  transfer 
parameters  h  and  7  are  bimodal  because  the  MCMC  eventually  finds  a  new  set  of  parameters  with  similar  error 
as  the  least-squares  optimal  parameters.  Further  investigations  are  necessary  to  determine  whether  relaxing  the 
bounds  is  appropriate.  For  control  algorithms,  we  are  mainly  concerned  with  bounds  on  the  parameters.  The 
95%  credible  intervals  for  all  the  parameters  are  listed  in  Table  3. 


5.  CONCLUSION 

In  this  paper,  we  employed  a  model  for  a  flexible  beam  actuated  by  a  single  SMA  tendon.  Additionally,  we 
investigated  a  systematic  approach  to  quantifying  parametric  uncertainty  using  a  Markov  Chain  Monte  Carlo 
algorithm.  The  results  indicate  that  the  posterior  densities  (a  measure  of  model  uncertainty)  appear  to  converge, 
but  additional  simulations  may  be  necessary  in  order  to  guarantee  convergence.  Additionally,  even  though  each 


a  [mm] 


Table  2:  Initial  optimization  values  and  results. 


Variable 

Lower 

bound 

Upper 

bound 

Initial 

estimate 

Optimal 

value 

Units 

Ea 

20.0 

90.0 

30.7 

20.0 

GPa 

Em 

20.0 

50.0 

26.0 

50.0 

GPa 

°L 

200 

350 

295 

350 

MPa 

/\(Jp 

7.0 

12.0 

9.2 

7.04 

MPa/K 

ep 

3.5 

6.0 

4.36 

5.11 

% 

T 

0.001 

2.0 

1.81 

1.2275 

s 

V 

lxlO-28 

lxlO-23 

l.Ox  10-25 

7.5xl0-24 

m 

h 

0.1 

3.0 

0.581 

1.418 

- 

7 

0.5 

5.0 

2.50 

4.62 

- 

£p 

1.0 

5.0 

4.0 

3.43 

% 

a 

0.1 

3.0 

1.0 

2.0 

mm 

L 

60 

120 

88 

106 

mm 

EI 

0.1 

5.0 

1.30 

0.99 

N-cm2 

(d)  (e) 

Figure  3:  MCMC  threads  (different  colors)  for  (a)  a,  (b)  sp,  (c)  h,  (d)  7,  and  (e)  EI. 


Table  3:  Comparison  between  95%  credible  intervals. 


Parameter  [units] 

a  [mm] 

eP  [%] 

h 

7 

EI  [N-cm2] 

Lower  Bound 

1.9746 

3.4006 

1.4310 

4.3674 

0.9123 

Upper  Bound 

2.0392 

3.4871 

1.6713 

4.9955 

0.9976 

(d)  (e) 

Figure  4:  MCMC  histograms  for  (a)  a,  (b)  Sp ,  (c)  h,  (d)  7,  and  (e)  EL 


parallel  thread  follows  a  different  path,  they  are  all  initialized  at  the  same  point.  It  may  advantageous  to  using 
different  initial  values. 

Future  work  will  compare  different  proposal  functions  and  their  effect  on  the  convergence  of  the  posterior 
density.  Additionally,  MCMC  results  will  be  compared  to  other  uncertainty  quantification  approaches  such  as 
bootstrapping.  Ultimately,  the  results  will  be  used  in  robust  control  algorithms. 
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